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Abstract. The multiconfiguration methods are widely used by quantum physicists and chemists for numerical 
approximation of the many electron Schrodinger equation. Recently, first mathematically rigorous results were 
obtained on the time-dependent models, e.g. short-in-time well-posedness in the Sobolev space H 2 for bounded 
interactions 1201 with initial data in H 2 , in the energy space for Coulomb interactions with initial data in the 
same space [251 [5] , as well as global wcll-poscdncss under a sufficient condition on the energy of the initial data 
[4] [5] . The present contribution extends the analysis by setting an L 2 theory for the MCTDHF for general 
interactions including the Coulomb case. This kind of results is also the theoretical foundation of ad-hoc 
methods used in numerical calculation when modification (" regularization" ) of the density matrix destroys the 
conservation of energy property, but keeps invariant the mass. 



1. Introduction 

The Multi- Configuration Timc-Dcpcndcnt Hartree-Fock (MCTDHF in short) system aims to describe the 
time-evolution of systems composed of fermions (e.g. electrons), that is, particles with half integer spin. More 
precisely, it represents a hierarchy of approximations of the linear TV particle (fermion) Schrodinger equation. 

d 

(1) i—^(t,x\,... ,xn) = H N <&(t,xi,.. .,xn), 

where ^ is the N-particle wavefunction which we normalize to one in L 2 (R 37V ;C) since \^(t, x\, . . . , x n )\ 2 is 
interpreted as the probability density of finding the i th particle in the position Xi at time t for all 1 < i < N. 
In order to account for the Pauli exclusion principle which feature the fermonic nature of the particles we are 
dealing with, an antisymmetry condition is imposed to the wavefunction "J i.e. 

*(i, xi, . . .,x N ) = e(r) *(x r (i), . . . ,x T (i)), 

for every permutation of {1, . . . ,7V} with signature e(r). The TV-body Hamiltonian of the system is then the 
self-adjoint operator Hm acting on the Hilbert space L 2 (R 3Ar ; C) given by 

(2) h n= {-\k,+U{x 1 )\+ Vdxi-Xjl). 

l<i<N ^ ' l<i<j<N 

The subscript Xi of A Xi means that the derivation is with respect to the space position of the i th particle and 
U{x) is the "external potential", created, e.g. by nuclei localized at fixed positions R m with charges z m > for 
1 < to < M . The last term of Hn is the interaction potential between the electrons, which is fundamentally 
given by the Couloumb interaction. Hence we have 

M 

(3) tf (a; ) = _ £ _=L_^ and V(\x-y\) = - 

\x - Rm\ \x- y\ 

We use the so-called atomic units for writing the Hamiltonian Hn and the N-particle Schrodinger equation; 
hence the Planck constant H, the mass of the electrons, the elementary charge and the factor j^- are all set 
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to one with eo denotes the dielectric constant of the vacuum. For the sake of simple notation we omit the spin 
variable, taking the spin into account does not add any mathematical difficulties. 

A particular advantage over the simple Time-Dependent Hartree-Fock (TDHF) method (see e.g. [23 | \13 \ l6l 15] 
is that MCTDHF can handle "correlation", an essential concept of many electron systems that vanishes (by 
definition) for TDHF. We refer the reader to [HI [19] for more details. 

MCTDHF methods are widely used for numerical calculation of the dynamics of few electron systems in 
quantum physics/chemistry (see e.g. [TJ [5J 131 DS] and references therein, also [7J |5J GO [T5] for the MCTDH case). 
From a theoretical point of view, excellent results were obtained for the stationary case [HI 03J [52] . However, 
the time-dependent case still poses serious open problems for the mathematical analysis. In fact, even the 
global-in-time existence of solutions without recourse to extra assumption as in [4] [5] is not proved yet. 

The MCTDHF system for the approximate wavefunction ^Mc(t, %i, %n ) is composed of A > N non- 
linear Schrodingcr-type evolution equations (for "the orbitals" , as a dynamic basis for an expansion by "Slater 
determinants") coupled with r := first order differential equations (for " the coefficients" C). The many 
particle wavefunction ^ can be well approximated by such linear combinations of Slater determinants, much 
better than by the simple TDHF method that corresponds to the limiting case K = N. In principle, for fixed 
iV, the MCTDHF equations yield a hierarchy of models with increasing accuracy with increasing A, in the sense 
that many particle wavefunction constructed from the solution of MCTDHF converges (in some sense) toward 
the exact solution with increasing A. However, especially in the time-dependent case, there is no proof for 
this seemingly "obvious" property. 

Let us now formulate the MCTDHF equations which are much more complicated to write down than the "usual 
NTS" like cubic NLS or " Schrodinger-Poisson" . For a short and readable introduction to the multiconfiguration 
time-dependent Hartree-Fock system, we refer the reader to [501 HS1 SI [S], or (2T[ Q31 [55] focussed on the 
stationary case. First of all, we introduce the set Fn.k of "cocfficents and orbitals" (C, $) 



(4) T> 



N,K 



|c = (C CT ) CTeEjViK eC, J2 l^l 2 = 1 | 



(5) |$ = (0i,..., 4>k) eilKT, J 0,0,(1.,- 

with 8i j being the Kronecker delta, the bar denotes the complex conjugation and 

H N , K = {a = {a(l) < . . . < a(N)} C {1, . . . , A'}}, \E NiK \ = := r. 

That is, the range of the family of increasing mappings a : {1, . . . , N} — > {1, . . . , A} for 1 < N < K with 
integers N and A and we shall not make difference between mappings and sets a for simplicity. Now, given 
a G S n.Ki we define the associated Slater determinant as follows 



^ a (xi,...,x N ) = 



4>a{N){xi) ■•• 4><j{N){xn) 

That is, the determinant built from the fa's such that i G a. The fa's being an orthonormal family, the factor 
^= ensures then the normalization ||$o-||z,2(R3iv) = 1. The MCTDHF wavefunction reads now 

(6) *Mc(t,x l ,...,x N ) :=tf M c(C(t), $(*))= 

,Xi, . . . ,X N ). 
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In the sequel, we shall use the same notation for operators acting on L 2 (R 3 ) and diagonal matrices of operators 
acting on vectors in L 2 (R 3 ) K . The equations of motion associated to the ansatz ([6]) correspond to a variational 
procedure combined with space rotations thanks to the gauge invariance of the system we are dealing with. We 
refer to 120] and particulary to [5] . In our formulation we obtain the following coupled system : 



(7) S 



if t C(t) = K[*](i) C(t), 

i T[C(t)} f t $(t, x) = T[C(t)] [-±A X + U{x)] $(t, x) + (I- P#) W[C, x) $(t, x). 



Let us now define the different operators and matrices involved in the system S. First, K[$] denotes an r x r 
Hermitian matrix depending only on the orbitals <pi and the binary interaction V and is indexed by the sets 
er, t <E ^n,k as follows 

(8) K[<&] CT , T (t) = 1 0- ~ - *r\{i,J},a\{fc 1 I}(-l)^ (-l)L X 

where the symbols (— = ±1 are not relevant for analysis and are given by 

(10) (-!)*! = < 

[ (_1)<T- I (fe)+<T- 1 (0 if > Z, 

with cr _1 (fc) being the position of the entry k in the set er. 

The matrix r[C(f)] is a K x J\ Hermitian "density matrix" depending only on the coefficients as follows 

(11) iW)k;= E (-ir" 1(i)+r " w ^(i)c T (t). 

<AW=t\{j} 

The projection operator P$ is given by 

K 

(12) p $ (V,) = V/ ^M^Mdx&fox), 

i=i ^ 

for all ^ € L 2 (E 3 ). That is the orthogonal projector onto the space spanned by the <^s. Finally, the K x K 
Hermitian matrix W[C, $](t, x) is given by the entries 

(13) W[C, $] tf (x) = 2 V 7jM / fc (y) F(|x - y\) My) dy, 
with 

(14) 7^ = ^(1-M(1-M E (-%(- 1 )*,I^^r. 

<r\{i,j}=T\{k,l] 

Observe that the potential y appears in the definition of W[C, <f>] and K[$], but that neither the kinetic energy 
operator nor the interaction U appear in and hence in the equation for the C(t). We refer the reader to 
[5] for more details on this formulation that corresponds to a carefully chosen gauge. Hence we have to deal 
with a strongly non-linear coupled system of r = (^) first order ODEs and K Schrodinger-type PDEs. 
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A preliminary result on the existence and uniqueness of solutions to the MCTDH system S (i.e. the adequate 
version of the multiconfiguration models for bosons) has been established in j2Q] under the drastic simplification 
of bounded and smooth interaction potential V (and {7 = 0), where H 2 regularity for a solution associated to 
initial data in the Sobolev space H 2 was shown to hold up to a certain positive time. Indeed, for such potentials 
it should be straightforward to get L 2 solutions, too. 



For the case of the Coulomb potential, some results have recently been obtained in [25J [5] with less regularity 
of initial data, namely in H . However, all these results are local-in-time in the sense that the existence, 
uniqueness and regularity persist only as long as the density operator associated to the system remains of full 
rank K. That is, the matrix F[C(i)] remains invertible. In case of a "loss of rank" at a certain time T* , 
the well-posedness holds only locally- in-time, until T*. However, in [5J H] it is shown how the global-in-timc 
existence can be assured under an assumption on the energy of the initial state ^Mc(C(t = 0), $(t = 0)). 

The purpose of this contribution is, essentially, to establish an L 2 existence and uniqueness result of solutions 
for the Cauchy problem associated to the MCTDHF equations, the system 5. For the binary interaction V and 
the external potential U we allow for a large class including the Coulomb potential : 

(15) U(\x\),V(\x\) £Wi(R 3 ,M)+V»(R 3 ,R), Pl ,p 2 > |. 

In order to make the proofs less technical, we shall take from now on U, V G L P (R 3 , R) with p > | omitting the 
decomposition in two parts. One can easily see that the same proof holds mutatis mutandis in the decomposed 
case, for the price of technically heavier and longer formulation of the estimates. One more simplification will 
be used in the sequel, we will usually write U, V G L P (R 3 ) with p > |, however, this has to be understood in 
the sense that U G L P (R 3 , R) and V G L q (R 3 , R) with different p and q satisfying p, q > 3/2. 

Observe that the assumptions above hold true in the Coulomb case as can be seen using the following 
cut-off 

= 1 for 0<r<l ( V 1 (\x\) = ^1 e L P (R 3 ) 

(16) x(r) { G [0, 1] for 1 < r < 2 and set < 
= for 2 < r V 2 (\x\) = 1 ~j^J 3: ^ X °(R 3 ) 



with 



l<P<3. 



Now, let us introduce 
(17) dF N . K = I (C, $) G F N , K ■ rank(r[C]) = K 



With these definitions, we arc able to recall a result from [251 E] that will be useful in the sequel. 

Theorem 1.1. Let (C°,$°) £ dF N ,K H H\R 3 ) K be an H 1 initial data, U, V G L d (R 3 ) with d > §. Then, 
there exists a time T* > such that the MCTDHF system S admits a unique maximal solution (C(i), in 

C 1 ([0,T*);C)' r x (c 1 ([0,T*);H- 1 (M 3 )) A 'nC°([0,T*); J ff 1 (K 3 )) K ) 

depending continuously on the initial data [C , $ ). Moreover, for all < t < T* , we have 

• (C(i),*(t)) G dT N ,K, 

• (^Mc(C(t)^(t))\H N \^ MC (C(t)^(t))) L2iR3N} = (*Mc(C ,$ )\H N \y M c(C ,$ )) L2{M3N) . 
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Finally, either T* = +00 or T* < +00 and \\F[C(t)} 1 \\ 3 ^ 2 dt = +00. In particular 

lim sup IITO)]- 1 ^ = +00. 

This Theorem then yields a local well-poscdness result to the MCTDHF system. However, observe that the 
result is possibly global since, for the time being, there is no indication that T* is necessarily finite [2"5I [5]. 
As noticed before, the special case N = K corresponds to a single Slater determinant ansatz. In particular, 
the set Tn,n coincides with 8Tn,n and becomes a smooth manifold actually. In other words, the matrix F[C] 
reduces to a globally invertible matrix since it becomes the N x N identity matrix. Therefore, Theorem 11.11 
ensures the global existence of a unique solution to the TDHF system, that is, T* = +00 and it improves then 
previous results obtained in [13] . Actually, the proof of the energy conservation presented in [5] is much more 
readable than the one there. 

Finally, recall that this result (and the one we shall present in the next section) is valid in the case of the 
MCTDH. The so-called Schrodinger-Poisson system (SPS), which coincides with the Hartree system in the 
special case of "Bose Einstein condensation" when Coulombic interaction is used, can be also obtained as a 
limiting case of the MCTDH and our result applies obviously to this model, too. 

2. Main result 

In [55], an L 2 -Theory to a TDHF type model is established. However, that model is a peculiar mixed state 
formulation. Usually, Hartree-Fock is characterized by a finite number of TV equations with occupation numbers 
equal to 1 which is, of course, a completely different paradigm than a mixed state for a one particle model. 
That is, the eigenvalues of F[C(t)] are one for all time t > 0. However, the result obtained in [51] obviously 
also holds for finite N (even if the author does not explicitly remark that) and can obviously be adapted to the 
usual TDHF setting. In our work this TDHF result is improved. 

Independently, Castclla established an L 2 theory of the mixed state Schrodinger-Poisson system [11] . More 
precisely, he studied a system of infinitely coupled Schrodinger equations with self-consistent Coulomb potential. 
The initial data needs only an L? bound, so the initial kinetic energy can possibly be infinite. Moreover, he 
obtained a blow-up (resp. decay) estimates for the solution as time goes to zero (resp. infinity). 

Our work is inspired by Castella's results and Strichartz techniques. Our result on MCTDHF applies also to 
the TDHF and the TDH "pure state" case. We prove the following 

Theorem 2.1. Let the potentials U, V G L d (R 3 ) with d > | in the sense of fl5)) and let (C°, $°) G BFn.k be 
an L 2 initial data with full rank (cf Q), l[17\ )). 

Then there exists a time T* > (possibly T* = +00 ) such that the MCTDHF system S Q) admits solutions 
(C(t),®(t)) satisfying 

• Ce C 1 ([0,T*),C) r and $ e C°([0, T*), L 2 (E. 3 )) K . 

Moreover, for all 2 < q < 6 

i) <P(t) G LW^r ([0,T},Li{R 3 )) K . 

ii) The solution (C(t), is unique in the class 

L°°([0, T], C) r x L°°([0, T], L 2 (M. 3 )) K n L^ 3 ^ ([0, T], L q (R 3 )) K , 

for all T < T*. 



iii) (C(t),<S>(t)) G dF N , K for all t G [0,T*). 
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Hence, this result establishes an L? theory for the solution ^mc{C{€), 3>(t)) ([6]) of the MCTDHF equations 
as long as the density matrix r[C(£)] remains of full rank. Note that our theorem yields a global L 2 theory 
for the TDHF and the TDH models (the " pure state" versions of the models studied in [UJ as mentioned 
above) . 

The paper is structured as follows. In the next section, we collect some well-known tools, like Strichartz 
estimates and the properties associated to the semigroup U(i) generated by i^A on L 2 (R 3 ). Moreover, we 
prove a local existence result using a standard contraction argument in an adequate space X^, q for a given reals 
p, q and a nonnegative time T' . In section 3, we prove that this local solution satisfies an a priori estimate 
which will be crucial in order to prolongate the solution beyond T' , Section 4 is dedicated to the proof of the 
main result 12.11 Finally, we make some comments and application of this Theorem. 



3. A FEW TECHNICAL LEMMATA 



First of all, let us specify the notation we will use throughout this paper and recall some well known tools. 
The real p' will be the conjugate of p with 1 < p < oo, that is, i + y = 1. By abuse of notation, we 
denote LP = L P (K 3 ,C) but also LP = L P (M. 3 ,M.) when there is no confusion. The same notation will be used 
for L P (M. 3 ,C) K and will be specified explicitly when necessary. The associated norms will be denoted || • \\ LP . 
The same conventions are adopted for the Sobolev spaces LL 1 . (U(i)) tS K is the group of isometries (ei tA )teR 
generated by |A on L 2 (R 3 ,C). Finally, k will be an auxiliary positive constant depending on N and K. Also 
"const." will denote generic constants depending on quantities that will be indicated explicitly when necessary. 
Next, for a given let for all 1 < i,j < K 

Y[$}ij(t,x) := f 4> i (t,y)V(\x-y\)~$ j (t,y)dy 

JR 3 

and 



D v [f,g}(t)= I f(t,x)V(\x-y\)g(t,y)dxdy. 

JR 3 xR 3 



From this point onward, we shall omit the dependence on t and x when the context is clear. Finally, for a given 
T > we denote L p T q = LP([0, T],L q ). Now, 

Definition 3.1. The pair of reals {p,q) is said to be admissible, we denote (p,q) G A, if and only if the 
following relation holds true. 

A = ( _) an d 2<q<6. 
dp 2 q 

Then, we arc able to recall 

Lemma 3.2. Let < t < T, then for all (a, b), (p, q) e A, (j) € L 2 and Lp G L?p' b , there exists p(a) and p(a,p) 
such that 

(18) ||U(t) 0|U, < p(p) U\\ L2 and I / U(f - s) V (s)ds\ Lr < p(p, a) \<p\ , tt , . 

Jo 

Proof. This is a classical Lemma and we refer the reader to [TBI LX3 [Ml HI] • D 

The first inequality appearing in (|18|) describes a notable smoothing effect. In particular it tell us that for 
all t £ K and (f> G L 2 , we have obviously TJ(t) <fi G LP . The second inequality is crucial when dealing with 
non-linearities in the framework of Schrodinger type equations. Indeed, without loss of generality we write the 
following generic Duhamel formula 



4it) = U(t)<f>-t f U(t-s)f(i>(s)), 
Jo 
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for a given functional /. Then, the first inequality of (|18[) allows to control the L 2 norm of U(t) (f> in terms of 
the L 2 norm of <f>. However, it is merely impossible to control the L 2 norm of J Q U(t — s) f(ijj(s)) ds in terms of 
the one of ip for general non-linearities /. The second inequality of (fl"8j) will, then, enable us to control the L^ q 
norm of J Q * U(i — s) f(ip(s)) ds for a given T > and a couple of reals (p, q) € .4 in terms of the ,<z norm of 
/ and allows to conclude. 

The Duhamel formula associated to the MCTDHF system S for a given initial data (C°,$°) is written as 
follows for all time t such that T[C(i)] is invertible, 



(19) 



~c(ty 




- c o - 


*(*). 




U(i) $° 




•/' 






-/a 


U(i-s 



[*(*)] C(s) 



ds. 



Remark 3.3. The potential U being time-independent, we chose for simplicity to add it to the non-linear part. 
An alternative way to proceed consists in considering the linear PDE i-^u(t, x) = — ^Au(t, x) + U (x) u(t, x) and 
find adequate real p for U G LP such that one can associate to this flow a propagator that satisfies Strichartz-type 
estimates (fl"8|) . We refer the reader to, e.g., [28] . 



Next, formal functional analysis calculation leads to 
Lemma 3.4. Let U, V <E L d , (p, q), (p i( %) e A for i = 1, . . . , 4 and T > 0. Then, 



(20) 
(21) 



(22) 



3^1,1^1 



|V[$]i, 



2 93|| 



3 3 



1 



£>v(0i(i) ^ 4 (*))dt 



IV 4 -J--2 

< T 2 Zjfc=1 



m £ < n i&l™. 

£^4 



Proof. The proof is nothing but a straightforward calculation based on the well-known Holder and Young 
inequalities in space and time and we refer to [29] for similar estimates. □ 



An immediate corollary of Lemma 13.41 is the following 
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Corollary 3.5. Let U,V G L d ,d > § and (p, q = -^) G A. Then, for all T > 0, we have 
(23) WQltyj <«rH |tf| £d ||$|| LP , 5 , 



(24) 



T 

K[$(i)] di 



(25) |w[c,$] <&| i?x < kt!-' [v| £ „ \c\% r mi~.> \\m Lr , 

(26) |p*w[c,*] m^- < «t!-5 |y| Ld [$|»o. ia ||$||^. 

Proof. The proof relies on the estimates (|2ip and (f2"2"j) of Lemma 13.41 The assertion (|2"3"|) is easy, in fact one sets 
(Pii 9i) = (Pi <z) € .4 and gets the result. Next, given a potential V G L d , the matrix K involves elements of type 
Dv(4>i4'j] 4>k<i>i) for which one use the estimate (|22j) by setting, for instance, (p3, q^) = (pi, qi) = (oo, 2) G A and 
(Pi> <7i) = (P2,Q2) = (p> ?) G 4. Also, the vector W[C, $] $ involves terms of type V[$]j. that can be handled 
using the estimate (|2"T|) . In fact, we set for instance (p2)?a) = (P3793) = (°°;2) G 4 and (^3,93) = (p, ?) G 4. 
Finally, the vector P$ W[C, $] $ involves terms of type Dv(4>i4>j,4>k4'i) 4*1- Observe that Dv{4>i(t>j,(j)p4>i) is a 
time-dependent scalar, thus, estimating the left hand side in an L v T ,q leads automatically to an L q x norm on (pi 
in the right hand side. For convenience, we estimate this term in L^ 2 in order to get an ||<^| L «>,2 and use the 
same choice as in (1241). □ 



Remark 3.6. From the one side, observe that for d > |, we have obviously 2 < q = < 6. Moreover, the 
estimates of Corollarv \3.5\ involve T a with power a > so that T a — > as T — > 0. Indeed, since 2 < q < 6, we 
have < - — A < 1. This observation will be crucial in the sequel. From the opposite side, assume V bounded, 
that is, d = 00, then, the estimates of the Corollary 1 3. 51 are valid with q = 2 and a, the power of T, equal to 1. 



Next, given T > and (p, q) G A, we define the spaces 

5TP.9 _ r°o,2 n tVA yP,Q _ rr*r v crP.g 

Z-Irp J_Jrp I I J_Jrp ^ (\,rp /\ Z-Irp , 

endowed with the norms 

\4>\ zv = H\\ L ^ + uw^, \C,$\\ xr = \c\cr + m z ™. 

A topology on Zj! 9 and X^ q being defined, we are able to prove the following 

Lemma 3.7. Let U,VeL d ,d>^,(p,q)eA such that q = Then, for all T > 0, we have 

(27) \\KMt)]C(t)-K[^(t)]C'(t)\\ L1{[Q T]) < const! tM ||(C, $) - (C, 

(28) |tf (*-*') 1 <«r|-5||^|| id ||$-$'|| z| ,, ai 

(29) ||W[C,$] W[C',$'] <const 2 rf-5 ||(C, $) - (C, $')IU|-*> 

(30) ||P* W[C, $] $ - P$, W[C", $'] $'1^1,3 < const 3 |(C, $) - (C, 
wi#i consti, const2 and const3 depending on 

N,K, \\U\ Ld , \\V\\ Ld , \\<S>\\ L ~,z, \\&\\ L ~,z, \Cl Cr , \\C'\ Cr , \m\ L ,, q and |$'|U,. 
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Proof. For a given V, observe that the difference K[$] — K[$'] involves terms of type Dv{<j>i (j>j, 4>k <i>i) 



) that we manage as follows 



\D v (fa cj> v 4> k ft) - DvM y 3 A' k <f>'i)\ < \D v ([<t>i - 4> v 4>k <t>i)\ + \D v {4>'i - <t>i)\ 

+ |zv(# - + |iv(# ? if ^ 0, - Wi])\ 

Thus, thanks to (f2"3|). we get 

|K[$] C-K[$'] C"|| il([0)T]) < |K[$]| il([0 T]) C"|| cr + |K[$] -K[a>']|| L1([0 T]) ||C"|| c „ 
< «Tf"i \V\\ Ld ||C- C'|| cr 

+ K rf-4 ||y|| it! ic'ic [Mi-..!*!*™ + l^-Hl^li-] l*-*'li-. a 
rf-i l^L^V [i*i»o f . a |* , | LS ,.« + ||*'lI^.a||*L s ,.«] [* 

The proof of the estimate (|27| follows then by setting, for instance 



const i 



K \\v\\ Ld \mi^ \m\i^ 
K\m L 4c'\\ C r 



K \\V\\ 



\w\\ L ~A\n% 



A\n L ^ + w\l-Am\ L ^ 



Proving the inequality (|28| is obviously a straightforward application of (|23|) . Now, one has 

||W[C,$]$-W[C',$']#'L W < \\W[C, $]($-$')! 

+ ||W[C-C',$]$'| L ^ V 
+ |W[C,$-* , ]* , | i?r ^. 

Observe that W[C, $] is quadratic in C and <£. Then, by (f25j). we obtain 

|w[c,$]<k-w[c",$'] $'1^,,, < k tI-4 |v| i£i ||qic,| ( f>ll^^ll$- ( f>'IL- 
+ «r!-* \v\ Li i\c\v + \c'\ c Aml-4n»Ac-c'\u. 



■t«-> \\v\\ L 4C\\ 



'II 2 

c 



+ «<£% 



Next, setting 



const2 



« Ml- [iqi^li^ + |c"| 2 r [|*| i? . a + i*^-.,] Ki^.. 

«Ml< [\C\cr + \C'\cr]ml^W\L£ ; *- 



for instance, finish the proof of (|29p . It remains to estimate the projection part in L^,' 2 . For that purpose, we 

estimate first [P$ — P$<] £ in L^' 2 for a given function e i^' 9 and all (p, g) G A This can be achieved 

thanks to Holder inequality in space and time as follows 



(31) 



K 

[P* - P*'] Clz,^ < 51 [ ll(^^fc)llii([ 0l T]) l^fc ~ ^fclli~- a + \\(€>4>k - 4>'k)\\Li([0,T]) Wk\ L ™. 

fc=l 

< K 



I* - *'h~-° + K IKIU*' II^'IIl- 2 II* - 
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Now, we have 

||P* W[C,$] $ -Pi,' W[C",$'] &\ L i2 < ||P$[W[C,$]$-W[C",$']$'|| L i,,2 

+ |[P$ ~P$/] W[C,$'] $'1^.2, 

< k |w[c, $] $ - w[C, $'] $'| LP ^ 

+ K |W[C", $'] $'| l p'. 5 ' [||*||iP.< + ll^'llz,|f- 2 ]I* - &lz*«- 

The estimate above is an application of (J31]), first with $' = 0, £ = W[C, $] $ - W[C", $' and second with 
£ = W[C", $'] Finally, Let g = ^ and recall that ^ G L d with d > §. Then following ([29]), we get the 
desired estimate by setting for instance 

const 3 = const 2 |$| L ,.«|$| i? ., a + «||F|| i(J ||C'||^||$ / || i? ,, a ||$'||i ? , 2 [||$|| ir + l^'H^], 

which achieves the proof. □ 

Now, let i?, T > be arbitrary reals to be fixed later on. Moreover let (C°,$°) G ^jv,jr and introduce the 
closed ball 

(32) A*'«(ii) = {(C, $) e X™ : ||C, $|^p,, < i?} . 

This defines a complete metric space equipped with the distance induced by the norm of Xj.' q . Finally, introduce 
the following mapping 



(33) 





'C(-) 




- c o " 






K[$( S )]C( S ) 




7Tco $0 : 








-/ 








(is, 




_$(•). 




U(-) $° 


Jo 


U(-- S ) 


C/$(s) +L[C(s),$(s)] 







with 

L[C, $] := rp]" 1 (7 - P*) W[C, $] $. 

This formulation is then well defined as long as the matrix r[C(f)] is invertible. From now on, we will consider 
initial data (C°,$°) € 8Tn,k, that is, the associated first order density operator is of full rank. Thus, the 
quadratic dependence of T[C] on the coefficients C a and the continuity of the MCTDHF's flow guarantees the 
propagation of this property up to a certain time T* > 0, at least infinitesimal but possibly infinite. That is, 
r[C(t)] is of rank K for all t € [0,T*), hence invertible and we refer the reader to [5] for more details on this 
point. 

Now, we claim the following 

Lemma 3.8. Let U,V G L d ,d> §, (p,q = -M^) G A and (C°,$°) G dF N , K with |(C°, $°)| C r xi 2 < P- Then, 
there exist a radius R > and a time T with < T < T* such that the mapping n is a strict contraction on 
XP' q (R). Moreover, given (C"°,$'°) G dF N , K with \C'°, &°\ C r xL 2 < (3 ,then 

(34) |(C7,$) - (C>,&)\\ xr < const |(C°,<f>°) - $'°)| c . xi2 , 

where (C, $) and (C",$') denote the fixed points of ir associated with (C°,$°) and (C"°, < f>'°) respectively. 

Proof. The proof is based on the Lemma T3.7I By abuse of notation, ttco^o (C, $) will be used instead of the 
vertical notation (|33l) . Moreover for notation's lightness, we set 

S((j>)(t) := f XJ(t-s)<j>(s)ds. 
Jo 
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Next, let (a, 6) € A, (C*, $), (C, $') G T > to be fixed later on and t G [0,T]. Finally we set 

(p=f,Q=&)€A. Then 



||7r c o i$ o(C(i),$(t)) - 7r c ,o i$ ™(C"(t),$'(,^ IICr><i a 



< 



(35) (l + p(a))|(C°,* )-(C rt, ,$' )|| Crxia 

+ | /' [K[$( S )] C( S ) - K[$'( S )] C(a)] ds\ cr + \S(U*-U&) (t)| L ..» 

+ |S (L[C, $] - L[C", $']) (t)!^,, := 7] + . . . + %. 

The term 71 is due to the Lemma I3~2"1 More precisely, the first assertion of (JTSJ). Next, thanks to the inequality 
(12711 of Lemma 13.71 we have 



(36) T 2 < const! T?-' ||(C, $) - (C, 

Now, we use the second assertion of (fT8]) in order to estimate 73 and 71- We start with 

% < p(a,p)\U$-U&\\ L ^, q ,, 

(37) < p(a,p)KTl-ilU\\ Li \\(C,$)-(C',&)\\ x ^. 

The second line above is due to (|28p where we upper-bounded obviously the Z^ q norm by the X^ q one. Finally 

T 4 < \\s (riq-^ic, *] $ - rfc'^wic", $') (t)I £ ..* 

+ |5(p*r[c]- 1 w[c,*] ^-p^rfCj-^c",*'] $') (t)^..*, 

< p(a, P ) irtGj-^c,*] s-ric'^wfc",*'] 

+ p(o,p) |p*r[c] _1 w[c,$] p^ric'p 1 w[C,$'] $'|| L ^. 

Next, observe the trivial algebraic relation 

\\ric}- 1 -tic'}- 1 ! < |r[c]- 1 (r[c']-r[c])r[C]" 1 |, 

< K\T[C\- 1 \\T[C']- 1 \{\C\ + \C'\)\C-C'\ Cr . 
Thus, by the mean of (|25l26l29l30p . we get 

(38) 7~ 4 < p(a,p) consUT«-i - 
More precisely 

const 4 = /.[Fll^lirici-^iiricr'lic'l^licic, + |c"| cr ][i + is'li-.'!*'!*™] x 

x l^l^..,!*'^,., + ir^-^tconsta + consta]. 
Summing (|36M38|) and adding the sum to the first line of (|31)|) leads to 

UcoMC®, *(*)) - tt W o(C" (t), $'(t))« Crxi M < (1 + p(a))|(C°, $°) - (C"°, $'°)| c , xi2 

(39) + [constx + p(a,p) const 4 + p(o,p) K \\U\\ Ld ] rf _ * |(C,$) - (C,* 7 )^.,. 

Now, the inequality p9[) holds for any admissible pair (a, 6) £ .A. Therefore, we write it in the special 
case (a, b) = (p,q) £ A and then in the case (a, b) = (oo,2) G A. Moreover, we set (C°,$°) = (C"°,$'°) 
and use \\C, $\ x *,<,, \\C, &\ x p,< < R since (C, $), (C, $') G i^ 9 (i?). Recall that T < T* since the ini- 
tial data is in OJ-n.k, that is, the matrix T[C] and T[C] are invertiblc thus there exists 8 > such that 
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iriCl-^o.p.T*). WiC'l'HL^^T*) < 9 - This Procedure leads to 

hc°,*°(C(t),$(t)) -7r c o,4o(C' (*),$' (t))]*p.. < const 5 ( K ,i?,0,p)Tf- 



(40) x I 

Equivalently, if we set (C",$') = (0,0) and use the fact that |C°, $ | Cl - xi2 < /?■ Again, after summation 

(41) \\Tr c »^o(C(t),$(t))\\ xr < const 6 ( K , R,6, j3,p) T~«~i . 
More precisely 

const 5 = 2 K R i \\V\\ Ld [5 + p(p)(26 2 R 4 + (26 + 7)8R 2 + 56)} + p(p) K \\U\\ Ld , 
const 6 = 2(1 + p(p)) + 2 K R 5 [\\V\\ Ld + p(p)6\\V\\ Ld (l + R 2 )} + 2Rp(p) K \\U\\ Ld . 

Thus, we choose R and T such that 

(42) T < 



1- 


R 


2q N 

6-9 






constg 


J 



That is 
(43) 



R consts — constg < 0. 



consteTa 2 < R and consts T a 2 < 1. 



Hence, by ([401 141j) , 7r c o $o is a strict contracting map on Xtp q . Eventually, it remains to prove the continuous 
dependence on the initial data (f34|) . Again, the essence is the inequality (|39|) . Let (C, $) and (C",4>') be the 
fixed points of ttqo^o and tt^/o $/o respectively. Let us write again the inequality for (a, b) = (p, q) and then in 
the case (a, 6) = (oo, 2) and sum up. One obtains 



IKco^o (C(t), $(t)) - TT^o^o (C'(t), $' 



\(C,*)-(C,$')\ X ,„ 



< 2(1 + p(p))|(C u , $ u ) - (C' u , $' u )II C .xl 2 + const 5 ( K , R,8,p) T*" 2 $) - (C, $' 
Finally, the estimate above with (|4"3"[) proves the continuous dependence on the initial data (|34p . 

Now, we propose the following crucial proposition 



ii. xi; 



□ 



Proposition 3.9. Let (C°, $°) £ dT N , K , U, V £ L d {R 3 ) with d > ^ and (p, g = ^-) G A T/ien f/iere exists a 
T(p, ||C|| Cr , |$| i2 , |$°| L 2) > and a umoue solution (C(t),$(t)) £ X^ q to 



(44) 



~c(ty 




- c o - 




*(*)_ 




U(t) $° 


Jo 



K[$(«)] C(s) 
U(t - s) \u $(*) + L[C(s), *(«)] 



In particular, the function part of the solution satisfies 

(45) <2p(p)||ci> || i2 . 

Proof. Thanks to Lemma [3.81 there exists a T" > for which the integral formulation (|44[) associated to the 
MCTDHF admits a unique solution (C(t), <£(*)) G A"|; g . The main point then is to characterize < T < T' for 
which the property ([4"5]) holds. The function part of the solution satisfies 

$(£) = U(t)*° - i 5(L[C, $])(*) 
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Let < r < T < T' where T is to be fixed later on and recall that |F[C(t)] 1 | i oo( r*) — ^- Next, we estimate 
the LP' q of the right hand side and get using the estimate (fTSj) of Lemma 13.21 



< pip) + \\s(u m)\z*" + tsinc]- 1 w[c,$] $)(t)| L? „ 

+ |5(T[C7(s)]- 1 P*W[C > *] *)(t)I £? .«, 

< p(p) + W$l L ^ + irtc]- 1 ! «w[c,$] 
+ irp]- 1 ! |p* w[c,$] $| L i,2, 

< i*°ii ia +««Tf-' \\u\ Ld m L ^+ K eTi-i mLd icf cr m 2 L ^m L ^ 
(46) + K eTi-^v\\ Ld \cf cr mi^m%, q . 

Next, we follow the argument of Tsutsumi in [27]. That is, we choose T so small so that there exists a positive 
number r] satisfying 



< 



f(r h T) :=p(p) m„-r, + K0 V T<-i (\U\ Ld + \\V\\ L 4C\\l [l + r, ||<%~ 

0<V<2p(p) \\$°\\ L 2- 
For that purpose, it is sufficient to choose T > such that 

(47) T < inf |(2«(9 [|C/| Ld + |^M<7& ^ft-,* [l + 2p(p) ||<i>V ] ) A , T'j . 
Then, let 

(48) 77o = miii{0<r ? <2p(p)|<i> | L2 ; /fa,T) = 0}. 
Now, let y(r) = ||$|| i? , a . Then, by @6J), we have for < r < T 

f y(r) < p(p) l#°l £a + *(?y(r) rf-i (|c/| id + |yMc& [1 + y(r) m L? ,* 

(49) V 

[ y(r = o) = o. 

Then, if T is chosen such that (gTJ holds, then by (gSJ |49j), we get 

y(r) <%<2p(p)|*°| La , for all re[0,T]. 
Thus, passing to the limit as r —> T, we get by Fatou's Lemma 

<2p(p) ||<I> || i2 , 

which is the desired inequality. □ 
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4. proof of Theorem 12.11 

This section is devoted to the proof of the main result Theorem 12.11 It relies mainly on proposition 13.91 and 
Theorem 11.11 which assures the existence and uniqueness of a solution, in short time, to the MCTDHF system 
with function part in the Sobolev space H . 

From the one side, let (C°, $°) G d^N,K by proposition 13. 9\ there exists a certain time T > such that the 
integral formulation (|44|) admits a unique solution (C(t), $(£)) in the space Af^'^. 

Now, one writes one more time the Duhamcl formula associated to the function part of the system 5, more 
precisely, the second component of (fT9|) . Next, one uses the estimates (|23|25|26|) of Corollary 13.51 . as in ([35]) 
with (C, $') = (0,0), and gets thanks to continuous embedding 

L«'(R 3 ) ^ ff-^R 3 ), for all (p, (?) G .A 

and the absolute continuity of the integral (the Duhamel formulation) that 

*eC([0,T o ],F- 1 ) thus $ g a,([0,T ],i 2 ) 

where the subscript „, stands for weak. 

From the opposite side, let (C°' n , $ -") G 0Tn,k H iJ 1 (R 3 ) if be a sequence of initial data for n G N such that 

(50) ||C ^,$°'"| CrxL2 < |C ,$°| Crxi2 , (C ' n ,$ a - n ) in C r xL 2 . 

The Theorem 11.11 guarantees the existence of a time T n > 0, the existence and uniqueness of a solution in the 
interval of time t G [0, T n ) such that 

(51) C n G C 1 ([0, T„), C) r , $"eC ([0,r„),ff 1 (R 3 )f. 

This solution satisfies the following conservation laws 

\\C n (t)l Cr = \\C > n \ Cr , l^(t)l L2 = \\^ n \ L2 , ViG[0,T„), 

for all ?i G N. Moreover, the solutions (C",$™) being continuously depending on the initial data, one can 
assume (up to the extraction of a subsequence) that the sequence T n converges towards T* where T* denotes 
the maximal time of wellposcdncss of the Cauchy problem associated to the system S with initial data (C°, $°). 

Then, by ([50| . (C°,<3>°) and (C°' n , 3>°' n ) are fixed points of the functions ttc° $° an d irco.n $o,n respectively 
on the same closed ball of for all (p,2 < q = < 6) G .A. Thus, one can pass to the limit as n — > +oo 
using the continuous dependence on the initial data ([34)) and (|50|) . We get 

(52) \C{t)\\ cr = ||C°| crj ||$(i)|| L2 = ||$0|| L2 , for all < t < inf {T , T*} . 

Hence, (J52J) with $ G C w ([0,inf {T , T*}),L 2 ) leads to $ G C([0,inf {T , T*}),i 2 ). Moreover, the fact that 
C G C^QO^nf {T , T*}),C) r follows from the continuity of t ^ K[$(f)] for all t G [0,inf {T , T*}) which is a 
consequence of the continuity of the ^s. 

Eventually, recall that the function part of this solution satisfies the a priori estimate (|45[) . In particular, the 
time To depends only on the constants N,K,p, ||C | C ,. and || < i )0 || i2 . Thus, one can reiterate the argument with 
initial data (C(inf {T , T*}),$(inf {T , T*})), (C(2 inf {T , T^}),$(2 inf {T , T*})), ... up to T*~. We now 
check the uniqueness of the solution to the integral formulation associated to the MCTDHF, namely (fT9|) . The 
uniqueness on [0, inf {To, T*}] is given for free by the contraction argument. Again, one can reiterate up to T*~ 
and get the uniqueness on the whole interval [0, T*). The points i) and ii) are straightforward and the last point 
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to make clear before finishing the proof of the Theorem 12. 1[ is the equivalence between the MCTDHF as an 
integral formulation and in the distributional sense. For that purpose, we refer, for instance, to (paragraph 
4: Proof of Theorem 2.2, more precisely, the uniqueness part of the proof). 



Recall that the Coulomb potential satisfies after cut-off A- G L d + L°° with d G [|, 3 [. Thus, for p 



2d 

d-l : 



have obviously (omitting the case d=^)3<q<6 which is the result obtained by Castella in [TTj for the SPS 
system. Moreover, notice that the estimates of the Lemma T3.7I hold true for bounded potentials. In particular 
in such case one can show that the mapping 



~c 













u<$> 



K[$] C 

TIC]' 1 (J-P*)W[C, $] $(s) 



ds, 



is locally lipschitz in C r x L 2 . This shows in particular the restriction of the results obtained in [20] since in 
this configuration we need only to prove the preservation of Tn,k by the MCTDHF flow as long as the density 
matrix is of full rank. 



5. SOME COMMENTS AND CONCLUSION 

For the time being, all the well-posedness results on the multiconfiguration models [2UJ HSJ LH 13 hold up to 
a certain time T*, possibly infinite, for which rank(F[C(T*)]) < K. In [H[5], we obtained with C. Bardos and 
I. Catto a sufficient condition that ensures the global well-posedness. More precisely, starting with initial data 
(C*°, $°) G F N ,K n fl^R 3 )* satisfying 

/*Mc(C ,*°)|ffjv|* M c(C ,$ )\ < min (*Mc(C,*)\ H *r\**ic(C,*)\ 

\ I L 2 (R3«) (C,<S>)eF N ,K-l \ I L 2( R 3N) 

then, the MCTDHF system S admits a unique global-in-timc solution (C(t), $(t)) in the class 

C 1 ([0, +oo); C) r x [C 1 ([0, +oo);H- 1 (M 3 )) K n C°([0, +oo); iJ^R 3 ))^) . 

However, such a condition is not possible for an L 2 theory since the energy of the initial data is not well-defined 
for (C°, ( I >0 ) G J~n,k- In order to remediate this problem in practice in numerical simulations, physicists and 
chemists use to perturb the density matrix [3 [10] (see also [20j). In order to illustrate their method we use, 
without loss of generality, the following perturbation 

T e [C]=r[C]+eI K , 

with I K being the K x K identity matrix. The MCTDHF system S © with F £ [C] instead of F[C] is then 
e— dependent. Then, observe that the L 2 theory we proved in the present paper is also global for S e since the 
time T of the Proposition 13.91 depends only on quantities that the flow S e conserves. Consequently, one can 
iterate the argument in time and cover the whole real line. Notice that S e does not conserve the total energy 
because of the perturbation. We are then able to claim the following corollary for 6> e 

Corollary 5.1. Let U, V G L d (R 3 ) with d > | and (C°,<!> ) G Fn,k be an initial data. Then, the perturbed 
MCTDHF system S e admits solutions (C £ (t), $ c (t)) satisfying 

• C e G C 1 ([0,+oo),C) r and $ e G C*°([0, +oo), L 2 (R 3 )) K . 

Moreover, for all 2 < q < 6 

i) * 6 (t)GL^([0,+oo),L«(R 3 )) x . 
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ii) The solution (C £ (i), is unique in the class 

L°°([0,+oo),C) r x L°°([0,+oo),i 2 (R 3 )) K r 

iii) (C £ (i), $ e (i)) G clF/v^ /or allt>0. 



iq 
3(5-2) 



([0,+^),L?(R 3 )) 



A global i? 1 theory on the system <5> £ is not possible since the total energy is not conserved, of course, since 
r[C] r e [C] _1 7^ Ik- We mention that energy conservation can be proved using a variational principle (cf [5]) 
or in a direct way by working only on the equations of the flow S ([26]). Moreover we do not believe that the 
perturbed energy 

TV 

E 

.i=i 



4 a, 

2 



U{ Xi 



is decaying during the time evolution. An alternative strategy will be to pass to the limit e — > 0, however this 
depends on wether or not the system S e satisfies a uniform (in e) estimate which is related to the non-decay of 
the perturbed energy. We hope to come back to this point in a forthcoming work. 
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